from pylab import *

d0 = loadtxt('../outputs/xe_Yp_0.00.dat')
a0 = d0[:,0]
z0 = d0[:,1]
xe0 = d0[:,2]

d1 = loadtxt('../outputs/xe_Yp_0.24.dat')
a1 = d1[:,0]
z1 = d1[:,1]
xe1 = d1[:,2]

fig = figure(figsize=(7,5))

plot(z0, xe0, 'b--', label='$Y_p=0$')
plot(z1, xe1, 'r--', label='$Y_p=0.24$')
ax = gca()
ax.set_xlim(z1.max(), z1.min())
ax.set_xlabel(r'$z$',fontsize=14)
ax.set_ylabel(r'$X_{e}$',fontsize=14)
ax.grid()

legend(loc='center left', fontsize=14)

tight_layout()

show()
